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Abstract 

A significant asymmetry in baryon/antibaryon yields in the central region of high 
energy collisions is observed when the initial state has non-zero baryon charge. This 
asymmetry is connected with the possibility of baryon charge diffusion in rapidity 
space. Such a diffusion should decrease the baryon charge in the fragmentation region 
and translate into the corresponding decrease of the multiplicity of leading baryons. As 
a result, a new mechanism for Feynman scaling violation in the fragmentation region 
is obtained. Another numerically more significant reason for the Feynman scaling 
violation comes from the fact that the average number of cutted Pomerons increases 
with initial energy. We present the quantitative predictions of the Quark-Gluon String 
Model (QGSM) for the Feynman scaling violation at LHG energies and at even higher 
energies that can be important for cosmic ray physics. 
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1 Introduction 



The problem of Feynman scaling violation has evident both theoretical and practical 
interest. In particular, this question is very important [H |2] for cosmic ray physics, 
where the difference from the primary radiation to the events registrated on the ground 
or mountain level is determined by the multiple interactions of the so-called leading 
particles (mainly baryons) in the atmosphere. 

Despite the lack of direct measurements of Feynman scaling violation for secondary 
baryon spectra in nucleon-nucleon collisions at energies higher than those of ISR, some 
experimental information from cosmic ray experiments seems to confirm [31 HJ [S] the 
presence of significant Feynman scaling violation effects. Now the LHCf Collaboration 
has started the search |6l |7] of Feynman scaling violation effects for the spectra of 
photons (tt^) and neutrons in the fragmentation region at LHC energies. 

In principle, the violation of Feynman scaling in the fragmentation region should 
exist due to the energy conservation, since the spectra of charged particles increase in 
the central region. However, no quantitative predictions can be made without some 
model of the particle production. 

The Additive Quark Model [H |9] predicts the violation of Feynman scaling in the 
fragmentation region due to the increase of the interaction cross sections. However, 
to make a description of the energy dependences of the spectra as a function of 
additional assumptions and parameters are needed. 

The QGSM [TUl E] allows the calculation of the spectra of secondaries at different 
initial energies in the whole xp region. The QGSM is based on Dual Topological 
Unitarization (DTU), Regge phenomenology, and nonperturbative notions of QCD. 
This model is successfully used for the description of multiple production processes 
in hadron-nucleon [121 [131 [IH [IS], hadron-nucleus jTHl [Ej, and nucleus- nucleus [IB] 
collisions. The quantitative predictions of the QGSM depend on several parameters 
which were fixed by comparison of the calculations to the experimental data obtained 
at fixed target energies. The first experimental data obtained at LHC [HI [20] show 
that the model predictions are in reasonable agreement with the data. 

We will consider the energy dependences of the spectra of secondary baryons in the 
projectile fragmentation region which we determine as the interval 0.05 < xp < 0.8. 
These values of xp are larger than the typical values for central production and smaller 
than the values where triple-Reggeon diagrams dominate. 

In the frame of QGSM several reasons for the Feynman scaling violation in the 
fragmentation region exist [21] • The first one is the increase of the average number 
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of exchanged Pomerons with the energy, which leads to the corresponding increase 
of the yields of hadron secondaries in the central region and to their decrease in the 
fragmentation region. This effect is present even at asymptotically high energies. The 
preliminary estimation of this effect was provided in [221 123] . 

In the case of nuclear (air) targets, the growth of the hN cross section with energy 
leads to the increase of the average number of fast hadron inelastic collisions inside the 
nucleus. Thus, the average number of Pomerons is additionally increased, resulting in 
a stronger Feynman scaling violation [2^ [25] . 

In [23] these predictions were taken into account to calculate the penetration of 
fast hadrons into the atmosphere, leading to a better description of the cosmic ray 
experimental data. 

The differences in the yields of baryons and antibaryons produced in the central 
(midrapidity) region of high energy pp interactions [13 [13 [2D1 [23 [23 [23 [2S] are 
significant. Evidently, the appearance of the positive baryon charge in the central region 
of pp collisions should be compensated by the decrease of the baryon multiplicities in the 
fragmentation region that leads to an additional reason for Feynman scaling violation. 
This effect has a preasymptotical behaviour and it is saturated at very high energies 
(see section 4). 

In the present paper we consider the effects of Feynman scaling violation, i.e. the 
energy dependences of the spectra of secondary protons, neutrons, and A produced in 
pp collisions in the fragmentation region. 

In our estimations the role of the nuclear factor for air nuclei should be similar to 
that presented in [221 123] . 

2 Inclusive spectra of secondary hadrons 
in the Quark-Gluon String Model 

The QGSM [TOl [TT] allows us to make quantitative predictions for different features of 
multiparticle production, in particular, for the inclusive spectra of different secondaries, 
both in the central and in fragmentation regions. In QGSM high energy hadron-nucleon 
collisions are considered as taking place via the exchange of one or several Pomerons, all 
elastic and inelastic processes resulting from cutting through or between Pomerons [21] . 

Each Pomeron corresponds to a cylindrical diagram (see Fig. la), and thus, when 
cutting one Pomeron, two showers of secondaries are produced as it is shown in Fig. lb. 
The inclusive spectrum of a secondary hadron h is then determined by the convolution 
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of the diquark, valence quark, and sea quark distributions, u{x,n), in the incident 
particles, with the fragmentation functions, G'^{z), of quarks and diquarks into the 
secondary hadron h. These distributions, as well as the fragmentation functions, are 
constructed by using the Reggeon counting rules [30j . Both the diquark and the quark 
distribution functions depend on the number n of cut Pomerons in the considered 
diagram. The details of the model are presented in references [101 [HI [121 [131 [U] • 




Figure 1: (a) Cylindrical diagram corresponding to the one-Pomeron exchange contribution to elastic 
pp scattering, and (b) the cut of this diagram which determines the contribution to the inelastic pp 
cross section. Quarks are shown by solid curves and the string junction by dashed curves. 

For a nucleon target, the inclusive rapidity (y), or Feynman-x (xp), spectrum of a 
secondary hadron h has the form [TU] : 

dn 1 da xe da ^ ,h/ a , j^h / \ 



dy ainei dy ai^ei dxp 



n=l 



where xp = 2p\\/^ is the Feynman variable, and xp = '^Ej^, and the functions 
determine the contribution of the diagram with n cut Pomerons and Wn is the 
relative weight of this diagram Y^^=iWn = 1- The last term in Eq. (1) accounts for 
the contribution of diffraction dissociation processes that are determined by the cuts 
between Pomerons {n = 0). 

For pp collisions 

(Ppp{x) = fgg{x+, n) ■ fq{x-,n) + /g Tl) ■ fgg{x-,n) + 

+ 2{n-l)f^{x+,n)-f^{x.,n), (2) 



Amj^/s + x"^ ±x 
4 



(3) 



where fgg, fg, and fs correspond to the contributions of diquarks, valence quarks, and 
sea quarks, respectively. 

These functions are determined by the convolution of the diquark and quark dis- 
tributions with the fragmentation functions, e.g. for the quark one can write: 

1 

fg{x+,n) = J Ug{xi,n) ■ G'^{x+/xi)dxi . (4) 

x+ 

The fragmentation functions G^{z) are independent on the number of cutted Pomerons 
n. On the contrary, the diquark and quark distributions u{x,n) (which are normalized 
to unity) become softer when n increases. Thus, for example in |T^ it was assumec^] 
that the diquarks distributions depend on n as: 

n,,(a;)~(l-a;)-"«+("-i). (5) 

If the intercept of the Pomeron trajectory is larger than unity, 

ap{0) = 1 + A , A > , (6) 

the average number of Pomerons which should be accounted for increases with energy. 
The probabilities for cutting different numbers of Pomerons, n, Wn in Eq. (1), can be 
calculated in the quasieikonal approach [31]. The results of the calculation at four 
different energies, -y/i = 17.3 GeV, 200 GeV, 8 TeV (the current LHC energy), and 
lOOTeV (significant energy for the Pierre Auger cosmic ray observatory, see, for example 
\32\ ) are presented in Fig. 2. 



3 Baryon/antibaryon asymmetry in the QGSM 

In the string models, baryons are considered as configurations consisting of three con- 
nected strings (related to three valence quarks) called string junction (SJ) [331 Ell EH 
136] . Such a baryon structure is supported by lattice calculations [37]. In the case of 
inclusive reactions the baryon number transfer to large rapidity distances in hadron- 
nucleon reactions can be explained [T51 EH ESI EHl SO] by S J diffusion. 

The production of a baryon-antibaryon pair in the central region occurs via SJ-SJ 
(SJ has upper color indices whereas SJ has lower indices) pair production, which then 
combines with sea quarks and sea antiquarks into a BB pair [351 SI] ) cis it is shown 
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Figure 2: The calculated probabilities for cutting different number of Pomerons at energies y/s ~ 17.3 
GeV (dotted curve), 200 GeV (dash-dotted curve), 8 TeV (dashed curve), and 100 TcV (solid curve). 

in Fig. 3a, the contributions of these processes to the inclusive spectra of secondary 
baryons being determined by Eq. (2). 

In the processes with incident baryons, e.g. in pp coUisions, another possibihty to 
produce a secondary baryon in the central region exists. This possibility is the diffusion 
in rapidity space of any SJ existing in the initial state and it can lead to significant 
differences in the yields of baryons and antibaryons in the midrapidity region even at 
high energies [15]. The most important experimental fact in favour of this process 
is the rather large asymmetry in fl and f2 baryon production in high energy 7r~p 
interactions |32]. 

The theoretical quantitative description of the baryon number transfer via SJ mech- 
anism was suggested in the 90's and used to predict [43j the p/p asymmetry at HERA 
energies. 

In order to obtain the net baryon charge we consider, following ref. [I5] three different 
possibilities. The first one is the fragmentation of the diquark giving rise to a leading 

^ There is some freedom [IQ in how to account for this effect. 
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Fi gure 3: QGSM diagrams describing secondary baryon B production by diquark d. (a) Central 
production of BB pair. Single B production in the processes of diquark fragmentation: (b) initial SJ 
together with two valence quarks and one sea quark, (c) initial SJ together with one valence quark 
and two sea quarks, and (d) initial SJ together with three sea quarks. Quarks are shown by solid 
curves and SJ by dashed curves. 



baryon (Fig. 3b). A second possibility is to produce a leading meson in the first break- 
up of the string and a baryon in a subsequent break-up (Fig. 3c). In these two first 
cases the baryon number transfer is possible only for short distances in rapidity. In the 
third case, shown in Fig. 3d, both initial valence quarks recombine with sea antiquarks 
into mesons M while a secondary baryon is formed by the SJ together with three sea 
quarks. 

The fragmentation functions for the secondary baryon B production corresponding 
to the three processes shown in Figs. 3b, 3c, and 3d, can be written as follows |15j : 

fliv ■ < ■ ^'-^ , (7) 
aN-vf,-z^-il-z) , (8) 
aN-e- vf^ ■ z^-"^J ■ (1 - zf , (9) 

where qn is the normalization parameter, and w^, f^, f£ are the relative probabilities 
for different baryons production that can be found by simple quark combinatorics 
[m 05] ■ Their numerical values for different secondary baryons were presented in [28] . 

The first two processes shown in Figs. 3b and 3c, Eqs. (7) and (8), determine the 
spectra of leading baryons in the fragmentation region. The third contribution shown 
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in Fig. 3d, Eq. (9), is essential if the value of the intercept of the SJ exchange Regge- 
trajectory, a^j, is not too small. In QGSM the weight of this third contribution is 
determined by the coefficient e which fixes the small probability for such a baryon 
number transfer to occur. 

In the case of pp collisions the most sensitive ratios to the values of parameters asj 
and e are the ratios oi B/ B m the central region. If the initial energy is high enough, 
one can neglect the contributions of the processes of figs. 3b and 3c, so the B / B ratio 
is determined by the contributions of figs. 3a and 3d. 

The spectra of antibaryons are described by QGSM [121 ttSl US] with reasonable 
accuracy, so the values of the parameters which determine the contributions of Fig. 3d 
can be extracted from the experimental data. 

The energy dependence of p/p produced in pp collisions in midrapidity region {ijcm ~ 
0) is shown in Fig. 4. The theoretical curve has been normalized to the experimental 
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Figure 4: The QGSM description of the energy dependence of the p/p in the midrapidity region. 
Solid curve correspond to the value asj — 0.5 and dashed curve to the value agj = 0.9. 

point at ^/s = 27.5 GeV, where the error bar is minimal and the analyses gives |27] : 

asj = 0.5 ± 0.1 , with e = 0.0757 . (10) 

Unfortunately, four experimental points at RHIC energy y/s = 200 GeV are in 
evident disagreement with each other, the calculation performed with values of the 
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parameters asj = 0.9 and e = 0.024 being in agreement with the lowest RHIC point 
(see dashed curve in Fig. 4). 

The data by the ALICE Collaboration JlBj for p/p ratios in pp collisions at -y/i = 900 
GeV and 7 TeV in midrapidity region are presented in Table 1, together with the 
QGSM predictions. These data confirm the asymmetry in the ratio p/p and they are 
in agreement with the QGSM predictions for the value asj = 0.5. 



SJ exchange 


^ = 900 GeV 


v/i = 7 TeV 


asj = 0.9 


0.89 


0.95 


asj = 0.5 


0.95 


0.99 


e = 


0.98 


1. 


ALICE 


0.957 


0.991 


Collaboration 


±0.006 ±0.014 


±0.005 ±0.014 



Table 1. The QGSM predictions for p/p in pp collisions at LHC energies and the corresponding data 
by the ALICE Collaboration. The value e = corresponds to the case without C-negative exchange. 

The LHCb Collaboration measured the ratios of A to A in the rapidity interval 
2 < y < 4 at i/s = 900 GeV and 7 TeV |17]. These preliminary data are also 
incompatible |18] with the QGSM calculation without SJ contribution {e = 0). Though 
the errorbars are too large to make any conclusion, the QGSM calculations with the 
value asj = 0.9 seem to be in a slightly better agreement with the data than the 
calculations with asj = 0.5. 

So, we can conclude that the transfer of baryon charge in large rapidity distances 
occurs up to the LHC energies, and it can probably be described by the QGSM by 
taking the values of the parameters asj and e presented in Eq. (10). 

4 Energy dependence of secondary proton spectra 
in the fragmentation region 

As it was mentioned above, the inclusive spectra of secondary net baryons B produced 
in the processes of Figs. 3b, 3c, and 3d are determined by the convolution of the diquark 
distribution u{x,n) in the incident particles with the fragmentation functions G^{z) 
presented in Eqs. (7) — (9). The baryon charge of all secondary particles is determined 
by the integral over these spectra and it should be excatly equal to two (i.e. to the 
baryon charge of the initial pp state). In the case the process in Fig. 3d is absent, the 
integral over the spectra of secondary net baryons totally saturates at the distance of 
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several units of rapidity from the projectile. Thus, in the absence of the process in 
Fig. 3d the Feynman scaling violation in the fragmentation region can only be due to 
effects connected with the increase of the average number of exchanged Pomerons. 

The SJ contribution to the inclusive cross section of secondary baryon production 
(Fig. 3d) at large rapidity distance Ay from the incident nucleon can be estimated as 

-^^as-e- e(i-"-)-^J^ , as = ■ vg . (11) 
(7 ay 

At asymptotically high energies, the baryon charge transferred to large rapidity dis- 
tances can be determined by integration of Eq. (11), and it turns out to be of the order 
of ^ 

{nsis oo))sj ~ as ■ jz r , (12) 

(1 - asj) 

only the left part of the initial baryon charge being available for the production of the 
leading baryons. 

The only free parameter in eqs. (7) — (9) is the normalization qn, that should be 
modified at asymptotically high energies in the presence of the SJ mechanism for the 
baryon charge as 

aN{s oo) = aN ■ -, — ^ —, — 7 TT — , (13) 

to guarantee the conservation of the baryon charge. One can see that {nB{s — > oo)) sj 
is finite if asj < 1 (see Eq (12)), so the Feynman scaling violation due to the discussed 
effects has a preasymptotical behaviour. 

To obtain the QGSM predictions for the spectra of leading baryons at finite energy 
s we have to calculate the value of (^b(s))sj at this energy for the renormalized 
value of aAr(s) in Eq. (13). This can be provided by the numerical integration of 
the convolution of the diquark distribution u{x,n) in the incident protons with the 
fragmentation functions G^{z). By this way we account for the rather complicate 
shape of (1/cr) ■ da^ /dy at small Ay. 

Though currently the value of asj = 0.5 seems more plausible [46j, the value of 
as J = 0.9 can not be excluded [171119]. Thus, in this paper we present the calculation 
obtained with these two values of asj, and also without any SJ contribution (e = 0.). 

The results of the calculations of the secondary proton spectra produced in pp colli- 
sions are presented in Fig. 5, together with experimental data [50j by the NA49 Collab- 
oration on proton spectra in pp collisions at 158 GeV/c beam momentum {\fls) = 17.3 
Gev) and the data [51j at 100 GeV/c (v^=14 GeV) and 175 GeV/c (v^=19 GeV). 
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Figure 5: The QGSM predictions for the spectra of secondary protons produced in pp collisions at 
-y/s = 17 GeV (left) and 8 TeV (right), compared to the corresponding experimental data by the NA49 
collaboration [50] at ^/s = 17 GeV (points) and data [51 at ^/s — 14 (triangles) and 19 (squares) 
GeV. Calculations without SJ contribution are shown by solid curves, results for asj = 0.5 by dashed 
curves, and for asj ~ 0.9 by dash-dotted curves. 



As it is seen on Fig. 5, the comparison of the QGSM calculations to the experi- 
mental data in ref. [50l IST] is rather good. The data of these two experimental groups 
are in good agreement at xp < 0,6 and in some disagreement at larger xp. The calcu- 
lated spectra of secondary protons produced at y/s=17 GeV and 8 TeV are significantly 
different. The calculated spectrum at ^/s = 17 GeV increases with xp rather monoton- 
ically, in agreement with the existing experimental data. This spectrum only weakly 
depends on the SJ contributions, the reason being that the value of {nB{s))sj at this 
energy is rather small. With the increase of the energy until y/s = 8 TeV the rapidity 
region accessible for the baryon number transfer increases, increasing also the value of 
{nB{s))sj and making the SJ effects more visible. 

The peak appearing at y/s = 8 TeV and xp < 0.06 is connected with the protons 
produced together with antiprotons via the mechanism shown in Fig. 3a. 

The proton spectra shown in Fig. 5 at both (low and high) energies depend on the 
values of several QGSM parameters. These considered values were determined from 
the comparison of the model calculations to the experimental data, but since there 
is some small disagreement between the calculations and the data. The main part of 
the uncertainty being connected with the normalization), the values of parameters, as 
well as the absolute values of the calculated spectra, can be known on the accuracy of 
the order of, say, 20%. Consequently, the difference between the calculations with and 
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without SJ contribution at the same energy can be estimated on the level of the 20% 
of this difference, in agreement with Eq. (10). 

Similar situation appears when we consider the ratios of the spectra of the same 
secondary h at different energies: 
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(14) 



In Fig. 6 the QGSM predictions for the ratio i?p(8TeV/17GeV) are presented. Here 
again the normalization uncertainties are canceled and the accuracy of the model pre- 
dictions is good enough. 
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Figure 6: The QGSM predictions for the ratios of the spectra of secondary protons produced in pp 
collisions at ^/s = 17 GeV and 8 TeV. Calculations without SJ contribution are shown by solid curves, 
results for asj = 0.5 by dashed curves, and for asj = 0.9 by dash-dotted curves. 



When the initial energy increases the spectra of secondary protons with xp > 0.3 
decrease, mostly due to the increase of the average number of cut Pomerons (see Fig. 2). 
The same effect increases the spectra at xp < 0.3. Some additional contribution to 
these Feynman scaling violations comes from the SJ effects which transfer the baryon 
charge to the Iow-xf region resulting in the renormalization of the parameter (see 
Eq. (13)). 
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It is also interesting to consider the spectra of secondary protons at fixed points 
of xp as functions of initial energy. Here again, the absolute normalization of the 
curve contains some uncertainties, but their relative change with the energy and their 
dependence on the SJ effects can be calculated with reasonable accuracy. The results 
of such calculations for secondary protons with xp — 0.7, 0.5, 0.2, and 0.05 produced 
in pp collisions at different energies are shown in Fig. 7. 
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Figure 7: The QGSM predictions for the spectra of secondary protons as functions of the energy at 
fixed values of xf- The left top bold curves show the spectra of net protons, i.e. the values of the p — p 
differences. For the four panel, solid curves correspond to the calculation without SJ contribution, 
dashed curves to the value asj = 0.5, and dash-dotted curves to the value asj = 0.9. 

These spectra increase with energy at comparatively low xp = 0.05 and 0.2 (except 
for the calculation with asj — 0.9 in the last case), and they decrease with energy at 
Xp — 0.5 and 0.7. This change in the energy dependence is explained by the growth of 
the average number of exchanging Pomerons (see Fig. 2). 
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The important feature is that at low xp (about 0.05 and less) one should discrim- 
inate between the total yield of secondary protons and the yield of net protons, i.e. 
the values of the p — p differences, by comparing the upper and lower sets of curves in 
the left top panel of Fig. 7. The spectra of antibaryons are not affected by SJ effects, 
so after subtraction of p spectra the SJ effects are more visible in the spectra of net 
protons. 

5 Predictions for the spectra of neutrons and A in 
pp collisions 

The LHCf Collaboration plans to investigate the Feynman scaling violation for neutral 
secondaries in the fragmentation region. In this perspective, the QGSM predictions for 
the spectra of secondary neutrons and A-hyperons produced in pp collisions at y/s — 8 
TeV are presented in Fig. 8. 
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Figure 8: The QGSM predictions for the spectra of secondary neutrons and of A-hyperons produced 
in pp colhsions at — 8 TeV. Calculations without SJ contribution are shown by solid curves, results 
for as J = 0.5 by dashed curves, and for asj = 0.9 by dash-dotted curves. 

The spectra of neutrons are similar to the spectra of secondary protons presented 
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in Fig. 5 but the fragmentation maxima are not so stressed. The spectra of A-hyperons 
do not present such a maxima due to the additional suppression of fast strange par- 
ticle production with respect to the non-strange secondaries, what leads to a faster 
decrease of uu and ud fragmentation functions into A at large z in comparison with 
the fragmentation into secondary nucleon. 

In both neutron and A production cases the SJ effects lead to similar corrections 
of the spectra than for secondary protons. 

The ratios i?„(8TeV/17GeV) and i?A(8TeV/17GeV) of the spectra of neutrons and 
of A-hyperons produced in pp coUisions, defined by Eq. (14), are shown in Fig. 9. 
These ratios show the expected effects of the Feynman scaling violation in different 




Figure 9: The QGSM predictions for the ratios of the spectra of secondary neutrons (left) and 
A-hyperons (right) as functions oi xp. The sohd curves correspond to the calculation without SJ 
contribution, dashed curves to the value asj = 0.5, and dash-dotted curves to the value asj = 0.9. 

,xj?-regions. The ratios for secondary neutrons are similar to the corresponding ratios 
for secondary protons presented in Fig. 6. The ratios for secondary A are slightly 
different due to the difference in the fragmentation functions. 

The energy dependences of the spectra of neutrons and A-hyperons at xp — 0.05 
and 0.5 are presented in Fig. 10. For both secondary neutron and A the spectra increase 
with energy at xp — 0.05 and they decrease at xp = 0.5. The predicted spectra of 
net baryons have a more complicated energy dependence at Xp = 0.05, at Xp = 0.5 
practically coinciding with the total baryon spectra. 

Sometimes it is more suitable to compare the spectra in the rapidity variable. In 
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Figure 10: The QGSM predictions for the spectra of secondary neutrons (top) and A-hyperons (low) 
as the functions of energy at fixed values of xf- The bold curves (lower sets of curves) on the two left 
panels show the spectra of net neutrons and A-hyperons. For the four panel solid curves correspond 
to the calculation without SJ contribution, dashed curves to the value asj = 0.5, and dash-dotted 
curves to the value asj = 0.9. 



Fig. 11 we present the spectra of secondary protons, neutrons, and A at ^/s = 17 GeV 
and 8 TeV, as functions of the rapidity measured from the beam, ybeam — y, defined 
in the cm. frame. If Feynman scahng would be preserved, the curves in the left-hand 
region of Fig. 11 (i.e. in the beam fragmentation region) should coincide at different 
energies, the appearing differences between curves showing the violation of Feynman 
scaling. All curves are terminated at the point ycm — 0. 
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Figure 11: The QGSM predictions for the spectra of secondary protons, neutrons, and A-hyperons 
produced in pp coUisions at -/s = 8 TcV (sohd and dashed curves), and at ^/s = 17 GeV (dash-dotted 
and dotted curves with asj = 0.5. 

6 Conclusion 

We present the QGSM predictions for Feynman scaling violation in the spectra of 
leading baryons due both to the increase of the number of cutted Pomerons with the 
energy and to the baryon charge diffusion at large distances in the rapidity space. The 
experimental search of the spectra of leading baryons at LHC should allow to confirm 
or discard these effects. The possible violation of Feynman scaling in the fragmentation 
region at very high energies was discussed in [6j, based on Monte Carlo calculations. 
In the QGSM hese effects exist, but rule, not numerically large. 

We have neglected the possibility of interactions between Pomerons (the so-called 
enhancement diagrams), since our estimations [52] show that the inclusive density of 
secondaries produced in pp collisions at LHC energies is not large enough for these 
diagrams to be significant. 

Concerning the LHCf project we can note that the multiplicity of A in the frag- 
mentation region should be of the order of 10—15% the neutron multiplicity, and so it 
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should be accouted for. The production of other hyperons should be several orders of 
magnitude suppressed. 

We have also to note that our results are in reasonable agreement with the calcu- 
lations in ref. 

A more detailed version of part of the calculations we include in this paper was 
already published in ref. |21j . 
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